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Abstract. - We study the collective behaviour of a close packed monolayer of non-Brownian 

particles at a fluid-liquid interface. Such a particle raft forms a two-dimensional elastic solid 

1-^ ' and can support anisotropic stresses and strains, e.g. it buckles in uniaxial compression and 

t~| , cracks in tension. We characterise this solid in terms of a Young's modulus and Poisson ratio 

Q ' derived from simple theoretical considerations and show the validity of these estimates by using 

O I an experimental buckling assay to deduce the Young's modulus. 

(N 
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fvj ' Introduction. - Particle covered liquid interfaces are increasingly being exploited in a 

[~^ , wide variety of technological and medical applications [1]. Coating a liquid drop with a 

^O ' hydrophobic powder renders the drop non-wetting with the resulting liquid marbles free to 

^^ I roll on rigid surfaces or even float on water [2], a feature that has also proven useful as an 

adaptation to life on small scales [3]. In a similar vein, it is hoped that by encapsulating 

the active ingredients of drugs within a monolayer of colloidal particles (thereby forming 

j^ ' a colloidosome [4]) more medicines will soon be administered by inhalation, so improving 

their efficacy. In the latter case, the particle coated interface exists only in the preliminary 

I ' stages of manufacture, but in neither case have previous investigations been concerned with 

■ O ' understanding the properties of particle monolayers at liquid interfaces. In this article, we 

show that in fact such particle monolayers behave collectively like 2-dimensional elastic solids 

and further, we characterize the properties of this two-dimensional solid. 

A variety of solid-like behaviours are observed for monolayers of hydrophobic particles 
sprinkled densely onto an air-water interface for a wide range (2.5 fiia - 6 mm) of particle 
sizes. For example, such a monolayer buckles under sufficient static compressive loading (see 
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H ', figs. n(a) and (b)) demonstrating that it can support an anisotropic stress. This stress state 

- - '■ can only be supported by a material with a non-zero shear modulus, which is the signature of 

a solid. Once the compressive stress is removed, the monolayer returns rapidly (~ 0(0. l)s) 

to the undeformed state, ironing out the wrinkles formed by the buckling. This elasticity is 

also reminiscent of a solid and is in stark contrast to what is commonly observed in both dry 
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Fig. 1 - Typical wrinkling patterns observed experimentally: (a) Polysytrene particles with mean 
diameter 300 /xm and (b) Lycopodium particles with mean diameter 30 //m. (c) shows the crack in 
a monolayer of particles with diameter in the range 150 — 212 fim near equilibrium. Here the needle 
is used only to place surfactant on the surface; the gradient in interfacial tension causes a crack to 
open. The scale bar in each image corresponds to 5 mm. 



and wet granular systems, which also resist compression but do not, on the whole, return to 
the undeformed state once the compressive stress is removed. 

Furthermore, the solid will also fracture under relatively small tensions and is manifested 
by the spectacular branching fracture induced by the addition of a drop of surfactant as shown 
in fig. n(c). Here, the disparity in surface tension between the clean and contaminated areas 
is enough to open up a tensile crack that propagates with surfactant advected to the crack 
tip by Marangoni flows. 

These behaviours are curious since the monolayer has solid-like properties only in the 
presence of the liquid and is reduced to a cohesionless powder once this liquid is evaporated 
naturally, making this an unusual composite material. Capillary forces are responsible for the 
formation of the solid monolayer via the aggregation of particles trapped at the interface and 
give the monolayer cohesion under deformation. The particle monolayers we study are dis- 
tinct from other two-dimensional interfacial systems such as Langmuir monolayers [5] because 
the large size of the particles (diameters of up to 6 mm) makes them non-Brownian objects 
which interact solely via capillary attraction. In this respect they are similar to conventional 
bubble rafts [6], and so we dub them particle rafts. However, bubble rafts of a similar scale 
are more easily deformable in compression and stack into multiple layers under the effects of 
large compressive stresses rather than forming larger scale wrinkles. This latter difference is 
due to the fact that bubbles are considerably less dense than the bulk fluid, ensuring that it 
is energetically cheaper to stack them than to deform the bulk interface significantly, which is 
a requirement of wrinkle formation. For a monolayer particle raft loaded in uniaxial compres- 
sion, the steric interaction between particles allows for anisotropic stresses while attractive 
capillary forces prevent stress relaxation due to particle creep in the direction orthogonal to 
the uniaxial compression. 

Although the dynamics of crack propagation in this system are of interest, we defer further 
discussion to the future and focus in this article on characterising the elastic properties of the 
particle raft. In this system, it seems reasonable to assume that the particle raft behaves 
approximately as an isotropic solid, in which case its elastic behaviour can be completely 
characterised in terms of the Young's modulus, E, and Poisson ratio, ;/; in the next section, 
we shall determine these parameters theoretically. 
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A Simple Model. - Microscopic theories of elastic solids usually rely on the existence of a 
mininuim in the interaction potential of the constituent particles, so that the elastic modulus 
is related to the curvature of this potential at the minimum. In the present situation, however, 
the following modification of a simple argument by Lucassen [7] suffices. 

For a conventional two-dimensional elastic solid, the mean stress, a = ^^^j^, and strain, 
e = ^^^j^ (where the subscripts 1, 2 denote the two principal directions), are related by [8]: 



2 



E 



(1) 



where v is the Poisson ratio of the solid. If we replace the stress by a thickness-averaged 
isotropic tension, r = ud, then this relation reads: 



We then have: 

l-v_de_ldA_ 1 d{Ai + A^) 

Ed "d^ ^ A'd^ ^ Ai + As d^ ^' 

where Ai is the area of the system covered by liquid. As is the area covered by solid particles 
and A = Ai + As is the total area of the system. The area covered by particles is constant 
regardless of the tension applied since the solid particles are effectively rigid at the individual 
level in comparison to the raft as a whole which collectively can be soft. In addition, Ai -jj- oc 7 
(with a constant of proportionality that may, in general, be different to 1) so that Q can be 
rearranged to give the effective Young's modulus: 

where (/> — Ag/A is the solid fraction of the interface. 

While naive dimensional analysis would predict that E ~ 7/rf, by recalling J^l, we might 
expect that E should only enter through the combination E/{1 — v). This would lead us 
to expect that i5 ~ (1 — v)f[(j))^/d for some function /(0), which we know must have the 
feature that as -^^ 1, / ^ cxd, since the particles are assumed rigid. The argument leading 
to Q showed that /(0) oc (1 — (?!))~^, which does have the desired divergence properties and 
is perhaps the most natural such function. 

An estimate of the value of v may be determined from the assumption that the particles are 
hexagonally close packed discs that interact only by means of an attractive central force and 
a repulsive steric interaction: a reasonable approximation for the case of spherical particles. 
Figure |21 demonstrates the origin of the Poisson effect for a single rhombic cell in the lattice 
in which the two central particles are displaced by a distance 2e. From elementary geometry, 
the rhombus joining the centres of the four particles in the undeformed state has width 2^/2>R 
and height 2i?, where R is the particle radius. During the deformation, the width decreases to 
2^R^/T^^2eJ^ « 2\/3i?(l - e/3i?) and the height increases to 2(e + R) so that the Poisson 
ratio is simply v = l/v^. 

Under the same assumption that the particles have a perfectly circular cross-section and are 
monodisperse, we may determine from elementary geometry that in hexagonal close packing, 
(j) — 7r/2-\/3. This gives an estimate of the typical value of the Young's modulus: 

E^2^P^1^Am1 (5) 
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Fig. 2 - Diagram showing the Poisson effect for four hexagonally close packed spheres (forming a 
rhombic cell) interacting via a purely central force. 



The calculations above are based on the assumption that the particles are hexagonally 
close packed discs when, in fact, their shape is irregular leading to non-ordered packings. 
However, this approximation provides estimates for the values of v and 0, which can then be 
used to obtain an estimate for the prefactor in lO 

The results (|4I5|) show that the elastic properties of the system are determined just by 
the surface tension, 7, and the particle diameter, d. In particular, there is no dependence 
on the particle elasticity since we have neglected the possibility of particle deformations in 
setting ^^ = in the argument leading to I0J. This is consistent with the picture that it 
is much easier to deform the layer by collective relative motion of particles than to deform 
the individual particles. We also see that the shape of the particles enters only via the 
dimcnsionless parameters 4> and v. 

Physically, the argument above shows that the elasticity of the particle raft arises because of 
changes in the interfacial area that are a result of the imposed deformation. In the undeformed 
state, the interfacial area is minimized (since this state is an equilibrium state) subject to the 
constraint that the solid particles are not deformable. Once a deformation occurs, the particles 
rearrange themselves to accommodate this imposition but the interfacial area associated with 
the new state is necessarily larger than that of the equilibrium state (since we have more 
constraints under deformation) and so energy is stored in the additional interface created. 

Experiment: a buckling assay. - To test the simple expression in Q, we used an elastic 
buckling assay of the particle raft to determine the wavelength (see fig. QJ of the instability 
and thence infer the value of the Young's modulus as a function of the particle size. To 
measure this dependence, each powder was in turn deposited by sprinkling onto an air-liquid 
interface contained in a home-made Langmuir trough (dimensions llcmx5cmxlcm for small 
particles sizes, 38.5cmxl8.5cmx4.2cmfor large particles). The barrier of the Langmuir trough 
consisted of a slide of the same width as the trough itself but slightly smaller thickness allowing 
the level of the liquid to equilibrate by flow beneath the barrier (see fig. EJ. The barrier was 
used to compress the particle monolayer until the buckling instability sets in; a schematic of 
the experimental setup is shown in fig. |21 The rate of compression leading to the buckling 
instability is not a factor in determining the wavelength of wrinkles produced as we measure 
the wavelength in the static configuration once any transients due to the induced flow in 
the liquid sublayer have decayed. Photographs of the wrinkled surface were taken using a 
high resolution CCD camera (1600x1200 pixels) for different runs. Single wavelengths were 
measured at different points of different photos (at least ten measurements) using image 
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Fig. 3 - Schematic of the experiment used to observe the buckling instability pictured in fig. 



analysis software (ImageJ, NIH) and the mean value calculated. 

The particles used were primarily made from Pliolite (Eliokem), a highly hydrophobic 
material, ground and sorted into different samples ranging in mean diameter from 40 //m to 
2mm. However, we also used a number of other hydrophobic powders both to obtain a wider 
size range and test the dependence of the results on the particles used. These included silanised 
Lycopodium [2] (grain size « 20 /im), graphite powder (grain size « 2.5 /im), powdered 
Chemigum (Eliokem, grain size k, 150 /im) and polystyrene powder of different mean grain 
sizes (ranging from 300 /xm to 7 mm) - the formation of a two-dimensional solid apparently 
only requires that these particles be trapped at the liquid-air interface. For particles larger 
than 40 /zm, we were able to use the CCD camera to measure the particle size and comparing 
20 grains of a particular sample gave us estimates for the polydispersity of that sample. For 
samples below this threshold, we relied on the size ranges specified by the manufacturer. 

Since the powders used were neither monodisperse nor uniformly spherical, the solid frac- 
tion of the monolayer, 0, was, in general, different from the maximum possible packing density 
of two-dimensional discs (which is close to 91%). Therefore, </> was measured just prior to the 
onset of instability using image analysis software (ImageJ, NIH) by changing the threshold 
on each image (the good contrast between liquid and particle phases ensuring that the exact 
value of the threshold used is immaterial). The experimental measurements of gave a mean 
value (\) = 0.85 with standard deviation 0.06. 

Most experiments were performed with the base liquid being water [p = 1000 kgm"'^ and 
7 = 72mNm~^) but we also performed experiments with saline solution (5%, by weight, 
NaCl so that p — 1034 kgm^'^ and 7 = 74mNm~-'^) to test for the possibility of electrostatic 
interactions mediated by the water. Previously it has been shown that the presence of the 
powder does not alter the surface tension coefficient as determined from comparisons between 
experiments on liquid marbles and numerical results for pure non-wetting drops [9] and so in 
the subsequent analysis we used surface tension coefficients for pure liquid samples. 

In fig. 01 we plot the wavelength as a function of the mean particle diameter for a range 
of different particle sizes and the two different base liquids. The results of these experiments 
lie on a single line and so we conclude that the precise geometry of the particles used as well 
as electrostatic effects are of minimal significance [10]. 

Determination of E from Buckling. - Classical elasticity theory [8] allows us to deduce 
the effective mechanical properties of this unusual two-dimensional solid from the measured 
wavelength which is assumed to be large compared to the particle size. We assume that the 
raft behaves as a thin sheet of isotropic, homogeneous elastic material of thickness equal to 
the particle diameter d in plane stress. At the onset of the buckling instability, the deflection 
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Fig. 4 - Experimental measurements of the dependence of the ripple wavelength A on the parti- 
cle diameter, d. All experimental points are for pliolite at an air-water interface unless otherwise 
specified by the key. The solid line shows the theoretical prediction 0. The error bars in the X- 
direction indicate the polydispersity of the sample, while the error bars in the Y-direction indicate 
the polydispersity in the measured wavelength. 

Fig. 5 - The Young's modulus as a function of the particle diameter for a range of particles. The 
solid line shows the theoretical relation E — 2.82"f/d. The dashed line shows the prediction from JSJ. 
All data points are for pliolite at an air-water interface unless indicated otherwise. 



of the raft is small and so its vertical displacement, h{x), measured relative to its equilibrium 
position in the absence of any compressive force satisfies the beam equation [8] : 






(6) 



where B = Ecfi /12{l — v"^) is the bending stiffness of the sheet, T is the compressive force 
per unit length in the sheet, p is the density of the liquid whose weight provides the restoring 
force, and g is the acceleration due to gravity. Substituting h{x) — j4sin(27ra;/A) into © gives 
us an expression for the compressive force T as a function of A; minimizing the result yields 
the selected wavelength A = 'k{4,EcP /{ipg{l — v'^))Y/^. Rewriting this expression gives the 
Young's modulus E in terms of the observed wavelength and the particle diameter as: 
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Using O the data shown in fig. ^ can be used to calculate the effective Young's modulus 
of the particle raft as a function of particle diameter d, with the result shown in fig. [S] We 
see that the Young's modulus is approximately inversely proportional to the particle diameter 
with a constant of proportionality close to the surface tension coefficient 7, as is predicted by 
our simple model 10}. 

Substituting the mean value for solid volume fraction measured experimentally (p = 0.85 
and V = l/\/3 into 10} we get E = 2.82j/d, indicated by the solid line drawn in fig. |S1 while 
the dashed line corresponds to |(S}- Linear regression performed on this data suggests that 
E = 0.154(i~^°^ in SI units, which compares favourably to the theoretical prediction based 
on = 0.85, which is -E = 0.205d~^. Finally, substituting the expression for the Young's 
modulus 10} with proportionality constant unity and rearranging |(2} yields an expression for 
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the wavelength of the wrinkles : 



1/4 
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where Lc = \J"i I pg is the capillary length, which result is plotted as the solid line in fig. 2|and 
fits the experimental data well. Linear regression on the data suggests that A — 0.235d°'*^, 
which is close to the theoretically predicted result, which in SI units reads A = 0.253d^/^. 

Conclusions. - We have seen that a two-dimensional monolayer of packed particles 
trapped at a liquid interface behaves like a two-dimensional elastic solid whose properties we 
have characterised. Our theoretical model for this behaviour shows that the properties of the 
particle raft are relatively insensitive to the details of the capillary interaction but do depend 
on the geometry of packing and has been verified by experiments on a wide range of particle 
rafts. Some variations in the experimental results could be accounted for by the polydispersity 
of the particle sizes and shapes and the resulting non-uniform packing. However, in such cases 
the value of v is different from l/vS, introducing an additional complexity that is beyond the 
scope of the present work. 

Many interesting questions associated with particle rafts involving both their dynamic and 
non-linear response remain to be answered. For instance, it is not clear what governs the speed 
of propagation of the cracks alluded to in the Introduction. Plastic behavior induced by shear 
in a Couette geometry shows that the particle rafts respond by localization of the deformation 
to a small region of the sheet, as for two-dimensional bubble rafts [11]. However our system 
may be more amenable to investigation since individual grains in particle rafts are effectively 
rigid. Indeed, this suggests that particle rafts could serve as a model two-dimensional granular 
system. 
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